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Abstract. The construction of an averaged theory of gravity based on Einstein's 
General Relativity is very difficult due to the non-linear nature of the gravitational 
field equations. This problem is further exacerbated by the difficulty in defining 
a mathematically precise covariant averaging procedure for tensor fields over 
differentiable manifolds. Together, these two ideas have been called the averaging 
problem for General Relativity. In the first part of the talk, an attempt to review 
some the various approaches to this problem will be given, highlighting strengths, 
weaknesses, and commonalities between them. In the second part of the talk, 
an argument will be made, that if one wishes to develop a well-defined averaging 
procedure, one may choose to parallel transport along geodesies with respect 
to the Levi-Cevita connection or, use the Weitzenbock connection and ensure 
the transportation is independent of path. The talk concludes with some open 
questions to generate further discussion. 



1. Smoothing of Space-time: What are the problems? 

The current widely accepted, standard cosmological model is a Universe that is 
homogeneous and isotropic on the largest of length scales filled with a mixture of 
baryonic matter, electromagnetic radiation, neutrinos, dark matter and dark energy. 
The detailed measurements of the cosmic microwave background via the WMAP 
program and other observational programs lend credence to this model. Interestingly 
enough, the two dark components are estimated to make up approximately 95% of 
the matter content of the Universe today and yet we have no true understanding 
of their composition or nature. Actually, the dark components can only be observed 
through their gravitational effects on baryonic matter and photons and are not directly 
observed. Can there be an alternative description for these observational effects that 
does not assume the existence of these mysterious dark quantities? 

One possibility is that both dark matter and dark energy are artefacts of some 
effective averaged theory of gravitation. If one employs Einstein's theory of General 
Relativity as our theory of gravity, then the resulting cosmological model based on 
the assumptions of homogeneity and isotropy is described by a single function of time. 
This idealized model is mathematically elegant and consistent with the observations 
on cosmological scales (provided one adds sufficient quantities of dark matter and dark 
energy to the model). However, can one simply ignore the structure and inhomogeneity 
that is present on smaller scales? Voids, walls, filaments, clusters, and super clusters 
are all examples of such structures that can be observed. Indeed, the smaller the scale, 
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the larger the inhomogeneity. Can one neglect the effects of these inhoniogeneities on 
our smoothed out idealized model? If not, how does one incorporate these effects into 
the model? Can the dynamical effects of these inhomogeneities be used as alternative 
explanation to dark matter and/or dark energy that is required in the standard model? 
Is General Relativity the correct theory of gravity to use on cosmological scales? 
The Einstein Field Equations (EFEs) of General Relativity 

Gapig) = uTais 

are currently used to describe the dynamics in both the realm of cosmology and 
the realm of the solar or stellar systems. However, most tests of General Relativity 
take place only on the scale of the solar system or stellar systems where Einstein's 
theory of General Relativity is consistent with all tests done to date 1 . Unfortunately, 
on much larger length scales. General Relativity appears to have some problems. 
Current gravitational and cosmological models using General Relativity fail to match 
observations on both the galactic scale (rotation curves of galaxies) and cosmological 
scales (acceleration of the universe) without the addition of the "dark" entities, matter 
and energy. 

What then might be the cause of the problem? What is different about the 
usual assumptions put into each model, solar system and cosmology. For solar system 
models, the right hand side is assumed to be an idealized, centralized, isolated source, 
however, for cosmological models, the right hand side is commonly modeled as a 
fluid. This means an averaging or smoothing procedure has been employed without a 
corresponding averaging or smoothing procedure on the left hand side of the equations. 
Shirokov and Fisher [2] had recognized the problem in the early 1960's, however it 
was not a well known issue until much later. It was only when George Ellis [3] gave 
a detailed description of the issues did the averaging problem receive the attention 
necessary to become a fundamental problem in mathematical cosmology. Both Ellis 
[3] and Shirokov and Fisher [2], suggested modified gravitational field equations for 
cosmology of the form 

Gap [g] = KTal3 = kT^p''^ 

where an overbar indicates that an averaging or smoothing process has taken place. 
However, this suggestion requires a response to the problem of averaging tensor fields 
on a manifold. Problem A: How does one Average tensor fields on a manifold? This 
is a non-trivial challenge due to the fact that there is no straight forward procedure for 
adding and comparing tensors located at different points on a differentiable manifold. 

Assuming that one has addressed problem A, one is left with an additional 
dilemma. Is it possible to relate G'q.^[5] with Gaplg]'^ Can we simply assume 
Gap[g] — Gq^[5]? In general, the answer is no, due to non-linearity of the EFEs. 
Both Shirokov and Fisher[2] and Ellis[3] introduce a Gravitational Correlation Tensor, 
Cap = Gai3 [g] — Gap [g] SO that the averaged field equations become 

Gap \g\ + Gap = nT^f^ 

but this leads us to one additional problem. Problem C: What is the nature of the 
gravitational Correlation , Cq^? 

Here we briefly review some of the approaches that have attempted to address 
problems A and C. 
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2. Proposed solutions to problems A and C 

2.1. Space-time Averaging 

In the 1960's, Shirokov and FisherP^ proposed a space-time averaging procedure of 
the form 

where x is the location of the macro-observer (center of averaging region) and ^ is 
the location of the micro-observer with respect to x. Unfortunately, the proposed 
averaging procedure is non-covariant, that is, the proposed average of a tensor is no 
longer a tensor. Further, in their analysis, the nature of the gravitational correlation, 
Ca/3 , was determined through perturbative techniques and therefore the gravitational 
correlation is non-exact. Consequently, the observations and results in this paper 
must be interpreted with care and consideration of these two weaknesses. 

The first covariant averaging or smoothing procedure appears to be found in a 
series of papers by Isaacson [4] in the late 1960's. Isaacson realized the problem of 
trying to add tensors at different points and therefore proposed an averaging procedure 
of the form 

T;{x)= I g2'{x,x')9l,{x,x')T^:ix')f{x,x')d'x' 

J all space 

where f{x,x') is a weighting function that falls to zero when x and x' differ more 
than some selected averaging scale and Jaii space d*^' ~ ^- quantity 

g^ {x,x') is the parallel propagator along the geodesic joining x and [3 [6]. The 
primary purpose of these papers was not explicitly for cosmology but to calculate the 
effective gravitational radiation in the high-frequency limit. Although the averaging 
procedure is covariant, the gravitational correlation was determined perturbatively 
and therefore the result is non-exact. It must be noted that in such an averaging 
scheme, the average of the metric is once again the same quantity, i.e, g^,^ = gajs. 
This observation leads us to question what geometrical object should be averaged to 
determine the effective averaged space-time. 

There have been a number of other space-time averaging procedures proposed 
since the 1980's. For example, Noonan [7] introduced the concept of micro and macro 
observers and the idea of duality. Unfortunately Noonan employed a non-covariant 
space-time averaging procedure and perturbatively determined the gravitational 
correlation. Therefore his results are both non-covariant and non-exact. Stoeger, 
Helmi, and Torres [5] employed a space-time averaging procedure which was equivalent 
to that of Shirokov and Fisher, hence the calculations are non-covariant. Since they 
applied the averaging procedure to a linearized model of gravity and a perturbed 
Friedman- Robertson- Walker model the results found therein are non-exact. 

In the early 1990's, Zalaletdinov [9] developed an averaging procedure that is 
similar to that proposed by Isaacson. To achieve a covariant averaging procedure, 
Zalaletdinov assumed the existence of bi- local transport operators A" {x,x') that map 
a tensor at x' to a tensor at a point a: in a manner that resembles the geodesic parallel 
transporter used by Isaacson. With some assumptions on the properties of the bi-local 
transport operators A" (a;, a;'), Zalaletdinov defined a covariant space-time averaging 
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procedure as 




l,^^^ A"^' {x, x')Al, {x, x')T-: (x') d^x' 



Instead of immediately averaging the EFEs, Zalaletdinov applied the averaging 
procedure to the Cartan structure equations. To determine the nature of the 
gravitational correlation, Zalaletdinov first defined a second order connection 
correlation tensor. With a few more assumptions, he was able to obtain splitting 
rules for various products of the Riemann tensor and the metric as well as a set 
of equations to determine parts of the connection correlation tensor. With these 
results, Zalaletdinov was able to apply the averaging procedure to EFE's to calculate 
the gravitational correlation in a non-perturbative manner. The resulting theory 
of Macroscopic Gravity is at present the only proposal in which the gravitational 
correlation is determined exactly and is based on a covariant averaging procedure. 
Particular cosmological models built on the theory of Macroscopic Gravity suggest 
that the gravitational correlation behaves as a spatial curvature term[Tni[II]. However, 
there are some weaknesses in the approach of Zalaletdinov. One of the weaknesses is 
the assumed existence and ill understood purpose of the bi-local transport operators. 
A second weakness results from the fact that the connection correlation tensor is not 
completely determined, and therefore there is still some ambiguity in the nature of 
the gravitational correlation jllj . The plethora of assumptions made to determine the 
gravitational correlation as found in the development of Macroscopic Gravity is also 
some cause for concern. 

2.2. Space Averaging 

Futamase [12] first studied the effects of averaging in a cosmological setting by 
employing a perturbative expansion of the metric and a simple spatial averaging 
procedure to determine the gravitational correlation. Futamase 's preliminary analysis 
is consequently non-covariant and non-exact. In a later paper, Futamase [T3] 
realizing the weakness in his previous work improved upon the averaging scheme by 
using the parallel propagator along a geodesic on 3-surface in a similar fashion to the 
work of Isaacson j^. While the procedure is not strictly covariant, the averaging of the 
3-geometry in this way is a covariant averaging process. The gravitational correlation 
was determined perturbatively, and hence the final result in this later work is also 
non-exact. However, even given the shortcomings in the overall analysis, Futamase 
determined that the dominant dynamical behavior of the gravitational correlation 
term behaves like a negative spatial curvature term in agreement with the conclusions 
found in [TUl [TT] based upon the covariant and exact procedure of Zalaletdinov . 

Kasai [14] approached the problem in a slightly different manner to that of 
Futamase |1 21 [T3] . Kasai foliated the space-time by flow orthogonal hyper-surfaces 
having 3-metric gij and concentrated his efforts on averaging a particular scalar 
quantity on the 3-dimensional hyper-surface. Kasai assumed an inhomogeneous and 
irrotational dust model with energy density p, and calculated the spatial average 
of the energy density which then defined the background energy density, pt, of the 
averaged space-time. With the assumption that the inhomogeneous dust model is 
an FRW dust model on average, Kasai was now able to define the scale factor a{t) 
of the background geometry via the relation SirGph — 3a^'^. With the background 
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geometrical quantities now known functions of the background energy density, Kasai 
defined two geometrical correlations [extrinsic curvature correlation and the Ricci 3- 
Curvature correlation] , and a density contrast function which completely characterize 
the gravitational correlation. Kasai determined the necessary conditions for the EFE's 
of inhomogeneous and irrotational dust models to have the form of a dust FRW on 
average. If these conditions were met, then some part of the gravitational correlation 
was determined. One of the strengths of the approach of Kasai is the fact that it is 
fully non-linear and exact. Unfortunately, the procedure applies only to irrotational 
dust models that are FRW on average. The weakness with this assumption is that the 
background geometry is fixed a priori and not determined via the averaging process. 
In addition, only part of gravitational correlation is explicitly calculated, some part is 
left undetermined. 

Buchert [HI [TB] improves upon the analysis of Kasai 14 by not fixing explicitly 
the background geometry a priori and concentrating his efforts on averaging scalar 
quantities only on the 3-dimensional hyper-surface. Buchert originally assumed 
that the cosmological model was an inhomogeneous dust model[15] foliated by fiow 
orthogonal hyper-surfaces with 3- metric gij , but later expanded the analysis to perfect 
fluid dust models [T6^. Buchert defined the spatial averaging operation 

T{X\t) = ^[ T{X\t)Jdet{g,,)(fX 
yo Jd * 

which is suitable for scalars on the 3-dimensional hyper-surface. Buchert defined 

the volume scale factor a^i (t) — \ ' 1 and applied the averaging procedure 
to the scalar parts of the EFE's. The two resulting equations now contain three 
unknown functions. The gravitational correlation is not completely determined. This 
observation is the primary weakness in the Buchert approach. Since the averaging 
procedure can only be applied to the scalar parts of the EFE's, the tensorial parts 
of the EFE's were ignored, and therefore only part of the gravitational correlation is 
determined in this way. However, given the weaknesses in the approach, the analysis 
indicates that the gravitational correlation term behaves as a curvature term [15] in 
some instances. 

2.3. Other Promising Approaches 

As an alternative to the approaches defined above, Boersma [17] defines a general 
averaging operator A and assumes that the Robertson- Walker geometry is a stable 
fixed point of this averaging operator. He further defines the linearized operator, A^^'^ 
for which stability of the averaging operation implies A°" [fSg^i/] — lim„_>.oo AS^''"'[5g^v\- 
He then shows that linearized averaging operation for metric perturbations, can 
be defined as a spatial averaging operation for scalars applied to Sg^Q and 5gl 
in synchronous coordinates. Unfortunately, the gravitational correlation is not 
completely determined. 

Employing the space-time averaging scheme of Zalaletdinov[9j and using a 
special limiting argument, Paranjape and Singh [18j successfully constructed a spatial 
averaging procedure that is based on a covariant and exact averaging procedure. 
Restricting their analysis to a cosmology in carefully selected gauges they observed 
that their equations agreed with, but were not the same as, the equations found in 
Buchert 's spatial averaging scheme for scalars described above. 
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More recently, Gasperini, Marozzi, and Veneziano 19_ proposed a Gauge invariant 
definition for averaging tensor fields. They introduce a Window Function similar in 
nature to Isaacsons weighting function f{x,x'). They make the argument that a 3 
dimensional spatial averaging can be determined from the 4 dimensional averaging 
scheme with an appropriate choice of the Window Function. Unfortunately, they did 
not apply the averaging procedure to EFE's and therefore have not attempted to 
address the nature of the gravitational correlation or problem C. 

In 2004, Debbasch [20] made arguments that an ensemble averaging procedure is 
necessary to average geometries. Unfortunately, no attempt to explicitly calculate the 
gravitational correlation resulting from averaging the EFEs is made. More recently, 
Sussman [21] defined quasi-local scalars in a class of inhomogeneous cosmological 
models, and was successful in averaging the scalar parts of the EFEs: in this way, 
it was similar in spirit to approach proposed by Buchert. Behrend [22) . has proposed 
a space-time averaging procedure based on the idea of a maximally smooth tetrad 
field which determines the averaged metric. Hehl and Mashoon [23]: developed a non- 
local theory of gravity based on a causal space-time averaged theory of gravity using 
GR|| which also shows some considerable promise. Khosravi, Mansouri and Kourkchi 
|24j have constructed some preliminary but important contributions to the ideas of 
"on" and "in" light cone averaging. Indeed, this idea has been further developed by 
Coley [3S] in which he discusses the need for Light-cone Averaging and averages the 
Raychaudhuri equation on the null cone in the same spirit of the 1-1-3 approach of 
Buchert. 

3. Another proposal to solve Problem A 

3.1. The Parallel Transport Problem 

Averaging involves integration/summation of tensor fields located at finitely separated 
points. This process is not straightforward on an arbitrary affinely connected and 
curved manifold. Specifically, one requires a notion of parallel transport of a tensor at 
a point x' along some curve Cx'~^x to a base point x to which we assign the average 
value in a unique and well defined manner. In general, the value of the tensor at x' 
when parallel transported to x is dependent upon the path Cx'^x, it takes to reach 
X, and the afRne connection used to facilitate the parallel transport. 

Two possible choices emerge that yield well defined transportation processes. The 
first choice is to parallel transport tensor fields along well defined curves with respect 
to a chosen connection, for example, parallel transport along geodesies with respect 
to the Levi-Cevita connection. Alternatively, and perhaps slightly more natural, one 
could select a connection such that the transportation process is independent of the 
path. Here we briefly investigate the two possible choices. 

Any tensor density field T with components T^^'"p" (x) can be parallel 
transported along a curve Cx'^x parameterized by z'^{t) from a point x' — z^{tq) 
to a point x = z^{t) via the elementary parallel propagators P^,, P^, and P^ 

rpai...a„f \_ p pQi pan pP'l pP'm rji^'i ■ ■ -a'^ / / \ 
^I3i...l3„, \^) — -^tu^a; ■ • ■ ^a',^ ^Pi • ■ • ^ft„ P'^-.-P'^, ^ > 

There is an excellent introduction to these elementary parallel propagators in the text 
by Plebanski and Krasinski [26] . Some of the properties of the elementary parallel 
propagators are described briefly below. 
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• The elementary parallel propagators satisfy the differential equations 



• The elementary parallel propagators have the fohowing initial values: 



• The elementary parallel propagators P^, P^, and P„ are functions of the 
connection, the endpoints x' and x and the curve C parameterized by z'^^t) 
connecting the two points. 

• Inverse Property: P;^,Pi' = S^, and P^,Pf = S^. 

3.2. Choice 1: Parallel Transport along Geodesic 

To begin, we must first select a unique curve that connects the points x and x' and a 
connection: for our purposes, we choose the geodesic and the Levi-Cevita connection. 
The geodesic is a "natural" choice as there are no other "natural" curves that connect 
x' and X. In Riemannian space, the geodesic is the shortest and straightest path 
connecting points x' and x. A weakness in this approach is the assumption that a 
unique geodesic exists connecting x' and x. The elementary parallel propagators 
no longer depend on an arbitrary curve and may be considered as functions of the 
endpoints x' and x. These special parallel propagators are usually denoted with a 
lower case g, i.e, gyj{x,x'),g'^,(x,x'),g'^ ix,x') to conform with standard notation for 
such objects O 

3.3. Choice 2: Path Independent Parallel Transport 

It is known that parallel transport is independent of path if and only if the curvature of 
the connection is zero. Therefore if one desires a path independent parallel transport, 
then one cannot use the Levi-Cevita connection to facilitate the parallel transport 
unless one is in flat space. One must employ a different connection, one in particular 
that has zero curvature. It is also known, that if the curvature is zero, then there exists 
n covariantly constant vector fields. Let e^" (i = 1, . . . ,n) he n linearly independent 
vector fields and assume they are covariantly constant with respect to some, as yet, 
undetermined connection. This requirement uniquely defines the affine connection 
W" = ei"e*/3.^ better known as a Weitzenbock connection. In this case, the 
elementary parallel propagators are factorable and can be shown to have the 
form 



PSir) =0 



p„(o) = i, p:,{o)^s:,, p:'{o)^s:' 





PS,{x,x') 



e^{x)e\,{x') 
e'„(x)e,"'(a;') 



5.^. Covariant Averaging Procedure 



Two options are proposed to define a covariant averaging procedure. The first is to 
parallel transport along the geodesic with respect to Levi-Cevita connection, in which 
case the curve and connection are now uniquely chosen. This was the approach used 
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by Isaacson^ in his investigations. The diSicuhy is that one is required to determine 
the geodesic first, and then exphcitly determine the parallel transporter via the initial 
value problem for the vector field to determine the nature of the parallel transporter. 
The second choice, is to choose an afhne connection that yields a Path Independent 
transportation process. Here, one parallel transports with respect to the Weitzenbock 
connection in which case the parallel transporter has a nice form and can be calculated 
quite easily. This approach closely resembles the approach used by Zalaletdinov in 
his investigations into averaging space-time [9]. Nonetheless, in either case, one can 
now integrate vector and/or tensor fields over compact regions of the manifold, and 
consequently, one can define an averaging or smoothing procedure. 

Definition 3.1 (Averaging/Smoothing Procedure) Let AA he a simply con- 
nected ^-dimensional manifold having a metric with Lorentzian signature. Let T^{x) 
be a continuous tensor field defined on some simply connected region TZ C Ai. Let T,^ 
be a compact subset of TZ at supporting point x. We define the average of the tensor 
field T'p{x), denoted Tp(x), as the definite integral at supporting point x, 



Note: one can substitute the geodesic transporter g'^,{x,x') if one decided to use 
geodesies and the Levi-Cevita connection. 

If one chooses to parallel transport with respect to the Levi-Cevita connection, 
then since the connection is metric compatible, we observe that 'g^ = gafj. That is the 
metric is an invariant under the averaging procedure. This observation may go against 
the natural idea that it is the metric that needs to be averaged. However, a rather 
nice observation results. The average of the Riemann curvature tensor corresponding 
to a space of constant curvature is also invariant under this averaging procedure. This 
implies that the average of a space-time of constant curvature, is itself, i.e., spaces 
of constant curvature are fixed points of the averaging/smoothing process, which is 
indeed a natural expectation of any averaging/smoothing procedure. If it does not 
make sense to average the metric to obtain a description of the averaged effective 
space-time, what then is the appropriate geometrical object of the microscopic space- 
time that should be averaged to determine an averaged effective space-time. Is it 
the Levi-Cevita connection? Possibly, but highly unlikely due to potential problems 
with the covariance requirement. Perhaps it is the curvature tensor Rp^g'^ A much 
better possibility? We have illustrated a covariant averaging procedure for tensor 
fields addressing problem A. We have not averaged the Einstein Field Equation's of 
General Relativity, and therefore have not addressed problem C of determining the 
gravitational correlation, so much more work to do. 

4. CONCLUDING REMARKS 

Averaging Spacetime: Where do we go from here? Before proceeding much further 
into this promising area of research, we, the research community, must address a few 
of the following issues. These are not in any particular order, but are all questions 
that have yet to be fully explored in a satisfactory manner. 




where 
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• On what length scale is Einstein's GR the appropriate theory of Gravitation? 

• If Einstein's GR is the appropriate theory describing the solar system and stellar 
systems, then what is the gravitational theory appropriate on cosmological scales? 

• Can a new theory of gravity be determined by simply averaging/smoothing the 
Einstein Field Equations? 

• Should the new smoothed theory appear as a totally new theory, or as GR plus 

new correction bits via a gravitational correlation? 

• Regarding the gravitational correlation (polarization) should it be determined 

— through perturbative techniques, or 

— assigned via geometrical assumptions a la Zalaletdinov? or 

— left undetermined as in Buchert? 

• Should an averaging scheme define an effective averaged geometry (as in 
Zalaletdinov) or should it define some averaged quantities defined on the un- 

averaged geometry (as in Buchert)? 

• Given that averaging/smoothing is much simpler to do in those spaces with tele- 
parallelism, should we average the GR|| equivalent to GR? 

• In cosmology, should wc ignore the effects of the micro geometry and concentrate 
solely on the macro, or should we attempt to merge both the micro and macro? 
For instance, in Maxwell's equations, the polarization terms are argued to be 
macroscopic in nature. 

• If the un-avcragcd space has a particular symmetry, then it makes some sense that 
any averaging procedure preserves these symmetries. How is this to be achieved 
in any averaging/ smoothing procedure? 

• One of our goals is to build better cosmological models. With this goal in mind, 
should we use a fully covariant space-time averaging procedure, or one better 
suited to cosmology, e.g., a 1+3 splitting of the space-time? 

• Can the inhomogeneities in the un-averaged geometry manifest an effective 
acceleration in the averaged geometry? This is the question that everyone would 
love to answer, but, this should come out of the model and not be a criteria that 
goes into the model. 

• Is the average of a bundle of null geodesies, a null geodesic? What about 
Causality? 

• Cosmology is tested with observations, and observations take place down the Null 
Cone: Should we not be averaging down the Null Cone? 

• Even if one is able to develop an averaging/smoothing procedure, and one 
is able to determine the gravitational correlation via some procedure and/or 
assumptions, one is still left with the Fitting Problem of How can one fit real-world 

observations to some smoothed out averaged model? 

Clearly, in order to proceed in developing a solution to the averaging problem, 
we must discuss and investigate a number of additional related issues concurrently. 
Hence the purpose of this session at the conference. 
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